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Abstract. The behavior of charmonia after the deconfinement transition is 
investigated on quenched lattices. Analysis of temporal correlators on fine lattices 
at temperatures upto 3 T c show that the J/tp and r\ c survive the deconfinement 
transition with little significant changes, and survive as bound states at least 
upto 2.25 T c . The spatially excited \ Ci states suffer serious system modifications, 
maybe dissolution, already a little above T c . 



The behavior of charmonia in heavy ion collisions has received considerable 
attention ever since Matsui & Satz suggested absence of J/ip bound states as a signal 
of deconfinement, based on nonrelativistic potential model arguments pQ. A different, 
more dynamical line of argument supported this suggestion, showing that the hard 
gluons available in the deconfined plasma will readily break the J/tp (treating J/tp 
as a Coulombic state) [2]. Direct studies on quenched lattices, however, show a quite 
different picture: the IS states like J/ip and rj c survive the deconfinement transition 
with little significant change, and bound states survive in gluonic plasma at least upto 
temperatures of 1.5 T c |310J[5]. Here we report an update of our results in reference 
PJ. Further details, including discussion of systematics, can be found in 

Lattice studies of finite temperature QCD employ the Matsubara formalism, in 
which field theory at thermal equilibrium is studied by trading the time axis for the 
inverse temperature. The detailed properties of mesonic states can be obtained from 
the Matsubara correlators 

G h (t,T) = (J h (t)4(0)) t . (1) 

Here Jh is the suitable hadronic operator, r denotes the Euclidean time and we 
consider zero momentum operators only. An integral equation then connects Gr (t, T) 
to the spectral function, crjj(o;, T), for the operator: 

r ( T\ f°w f ^ coshMT-l/2T)) 

Gh(t,T) = J o ^,T) sinh(cj/2r) ■ (2) 

The high temperatures that are of interest to us correspond to a small temporal 
extent. In order to have a reasonable number of data points for the correlator, we use 
very fine lattices. We present results here for two sets of lattices, with lattice spacings 
~ 0.02 and 0.04 fm, respectively, using the finer set to reach higher temperatures and 
the coarser set for temperatures closer to T c . For each set, we use a quark mass close 
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to the charm {mj/^ ~ 3.1 GeV for the coarser set and ~ 3.7 GeV for the finer set), 
and change the time extent N T to vary the temperature. Details of the lattices studied 
and the corresponding temperatures are tabulated below. 



a[fm] 


N% x N T 


T/T c 


# configs. 


a[fm] 


N'£ x N T 


T/T c 


# configs. 




48 3 x 24 


0.75 


100 




40 3 x 40 


0.9 


85 


0.04 


48 3 x 16 


1.12 


50 


0.02 


64 3 x 24 


1.5 


80 




48 3 x 12 


1.5 


60 




48 3 x 16 


2.25 


100 












48 3 x 12 


3.0 


90 



The ill-defined problem of inverting Equation J2J to extract a(uj,T) can be 
handled with the "maximum entropy method" (MEM) [7], which provides a prior 
guess for the solution through a Shannon- Jaynes "entropy" term(see 8 for a review) . 
The application of this method at finite temperature has led to important qualitative 
information about the dilepton yield from an equilibriated plasma 9 . However, the 
small extent of the Euclidean time direction at finite temperature makes the problem 
even more difficult 4 , and the role of the prior information becomes important (see 
for details). 

The pseudoscalar and vector states are explored by using the point-point 
operators Jh = cq$c and cj^c, respectively. The lightest states in these channels 
correspond to r/ c and J/ijj respectively. The extraction of the spectral functions below 
T c , using Equation @ and applying the MEM procedure, is stable and reproduces the 
properties of the ground state reasonably well. The spectral functions for these states 
at 0.9 T c (for the lattice with lattice spacing 0.02 fm; see table above) are shown in 
Figure ^ a )- Here we plot p(u),T) = o-(uj,T)/ui 2 . The ground state peaks and their 
strengths reproduce reasonably well the mass and amplitude obtained from a fit. The 
peaks at higher u> scale roughly with the lattice spacing and are probably dominated 
by lattice artefacts (see the discussion in reference 0). 
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Figure 1. a) The spectral function extracted from the temporal correlators at 0.9 
T c in pseudoscalar and vector channels. The first peaks correspond to rj c and J/tj} 
respectively, b) G(t,T)/G rccon,T* t 7 ", T), where G rcc0 n,T* (t - , T) is constructed 
from the spectral functions in a) (see text). 

As a first step towards looking for the effect of the deconfinement transition on the 
J I tp , one can ask the question whether the spectral functions in Figure ^ a ) explain 
the correlators above T c . We construct the model correlators above T c according to 
Equation (gj, 

cosh(w(r - 1/2T)) 



G 



rccon.T 



(r,T) 



du)<j(u), T*)- 



sinh(w/2T) 



(3) 
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where T* is a temperature below T c , in this case 0.9 T c . A comparison of 
G rec0 n,T* ( r j T") with the directly measured correlators G(r, T) give an indication of 
temperature modification of the mesonic properties above the transition. 

This simple comparison turns out to reveal a lot about the properties of J/ip and 
r\ c above the transition. Figure ^b) shows the result of such a comparison. For the r\ c 
channels, the spectral function at 0.9 T c is seen to completely explain the measured 
correlators at 1.5 T c , indicating that there is no significant change in this channel upto 
this temperature. At 2.25 T c , G(r, T) shows only small deviations from G roC on,T* (t, T), 
while at 3 T c significant changes are seen. For the vector channel, at 1.5 T c G(t,T) 
is described by G rccon ,T* (t, T) at small distances, while small deviations are seen at 
larger distances. The deviations appear at shorter distances at 2.25 T c , while pretty 
large deviations are seen at all distances at 3 T c . 

A more detailed view of the temperature modifications of the mesons can be 
obtained by extracting the spectral function directly at higher temperatures, by 
application of the MEM. At these temperatures, the small extent of the temporal 
direction makes a reliable extraction of the spectral function difficult without proper 
a-priori information. However, precise information about the structure of p(u>,T) at 
large ui, obtained from Figure^a), allows a reliable extraction of the spectral function 
UJ. As part of the prior guess, we provide the large uj structure from Figure ^a), 
smoothly connected to ~ uj 2 behavior at low u>. The spectral functions extracted with 
such a prior guess are shown in Figure [2J 
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Figure 2. Spectral functions in the deconfined plasma, for a) rj c and b) J/t/j 
channels. 

Figure [3 supports the trend seen in Figure Q]b). In the n c channel, a strong peak 
is seen at 1.5 T c , with an essentially unchanged position and strength from that below 
T c . In the J/ip also a strong peak is seen upto 1.5 T c ; no significant reduction in 
mass or peak strength is seen upto this temperature (see refernce 6 for a discussion 
of possible changes at this temperature). At 2.25 T c , a significant peak is still seen in 
both the channels, but with a much reduced strength. Furthermore, we do not see a 
statistically significant peak at 3 T c . (The error bars are standard deviations over the 
integrated strength, integrated over the marked region in the x direction; see reference 

While the IS states seem to be little affected by the deconfinement transition, the 
situation seems to be quite different for the IP states Xc and Xci ■ A similar analysis 
indicates that these states are significantly modified on crossing the transition point. 
In order to study these states, we use the point operators Jh{x) = cc and 0757^0, 
respectively. A reliable extraction of the spectral function below T c is possible, with 
the ground state peak reproducing the mass and strength obtained from a fit. At high 
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u>, lattice artefacts similar to those in Figure^a) dominate. In Figure[3]a) we show the 
comparison of the correlators for the axial vector and scalar channels at 1.1 and 1.5 T c 
with those reconstructed from the spectral function at 0.75 T c . The figure indicates a 
very significant modification of the properties of the states in these channels, already 
at 1.1 T c . Figure b) shows the spectral functions extracted above T c for the scalar 
channel. As above, we use the high energy structure below T c as part of the prior 
guess. The figure indicates that the \c state may have dissolved already at 1.1 T c . 
A very similar figure is also obtained for the \c% state. 
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Figure 3. a) G{t,T)/G rccon ( T 5 0.75T C ) for the scalar and axial vector channels, 
b) The spectral function for the scalar channel. 

These results have direct experimental relevance for the dilepton signal from an 
cquilibriated plasma. The differential dilepton rate is directly related to p(uj, T) in the 
vector channel: 

dW , 5a 2 1 , _ (AS 

Therefore our results for spectral function indicate that the dilepton signal should 
have a significant J/tp peak even for plasma temperatures above 2 T c . On crossing 
T c , the absence of the \ Ci states will cause a reduction of the J ftp peak due to the 
absence of the "indirect" J/ip from \ Ci decay (this may happen even below T c for 
full QCD). No further significant reduction of the peak strength will be seen at least 
upto temperatures of 1.5 T c , while reaching even higher temperatures will result in a 
weakening of the peak (possibly due to collision broadening). The results presented 
here, of course, only refer to a gluonic plasma. However, in a deconfined quark-gluon 
plasma the collision of the J/ip with the hot thermal gluons is expected to be the 
significant mechanism for dissolving the J/ip 2 ; the above results thus should hold 
at least qualitatively also for full QCD. 

This research is supported by BMBF under grant no. 06BI102. 

References 



Matsui T and Satz H 1986 Phys. Lett. B 178 416. 
Kharzeev D and Satz H 1994 Phys. Lett. B 334 155. 

Datta S et al 2003 Nucl. Phys. B (Proc. Suppl.) 1 19 487 (2002 Prepnnt |hep-lat/0208012) . 
Matsufuru H et al 2002 Preprint hep-lat/0211003 (Eur. Phys. J. C, to be published) 
Asakawa M and Hatsuda T 2004 Phys. Rev. Lett. 92 012001 (2003 Preprint hep-lat/0308034 I. 
Datta S et al 2003 Preprint hep-lat/0312037 (Phys. Rev. D, to be published). 
Asakawa M et al 2001 Prog. Part. Nucl. Phys. 46 459. 
Jarrel M and Gubernatis J E 1996 Phys. Rept. 269 133. 
Karsch F et al 2002 Phys. Lett. B 530 147. 



